OCU-PHYS-180 
January 2001 



Gauge-boson propagator in out of equilibrium 
quantum-field system and the Boltzmann equation 



A. NiegawaQ 

Department of Physics, Osaka City University 
Sumiyoshi-ku, Osaka 558-8585, JAPAN 



Abstract 

We construct from first principles a perturbative framework for studying nonequilib- 
rium quantum-field systems that include gauge bosons. The system of our concern is 
quasiuniform system near equilibrium or nonequilibrium quasistationary system. We 
employ the closed-time-path formalism and use the so-called gradient approximation. 
No further approximation is introduced. We construct a gauge-boson propagator, 
with which a well-defined perturbative framework is formulated. In the course of 
construction of the framework, we obtain the generalized Boltzmann equation (GBE) 
that describes the evolution of the number-density functions of gauge-bosonic quasi- 
particles. The framework allows us to compute the reaction rate for any process 
taking place in the system. Various processes, in turn, cause an evolution of the 
systems, which is described by the GBE. 
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1 Introduction 



Ultrarelativistic heavy-ion-collision experiments at the BNL Relativistic Heavy Ion 
Collider (RHIC) have begun and will soon start at the CERN Large Hadron Collider 
(LHC) in an anticipation of producing a quark-gluon plasma (QGP) [see, e.g., [[[], Q]]. 
The QGP to be produced is an expanding nonequilibrium system. Studies of the 
QGP as such have just begun. 

In previous papers, a perturbative framework has been formulated from first prin- 
ciples for dealing with out-of-equilibrium complex-scalar field system ||, O(N) linear- 
sigma system [|J, and the system that includes massless fermions |J. In this paper, 
we take up the out-of-equilibrium quantum-field theories that includes gauge bosons. 
We deal with the gauge-boson and the FP-ghost sectors in Coulomb gauge and, fol- 
lowing the same procedure as in [|, [|, §], we construct the gauge-boson and FP-ghost 
propagators and, thereby, frame a perturbation theory. Only approximation we em- 
ploy is the so-called gradient approximation (see below). We use the closed-time-path 
(CTP) formalism |], |7|, |J of nonequilibrium statistical quantum-field theory. It turns 
out that the form of gauge-boson propagator is quite complicated. 

Throughout this paper, we are interested in quasiuniform systems near equilibrium 
or nonequilibrium quasistationary systems. Such systems are characterized by two 
different spacetime scales |7|]; microscopic or quantum- field-theoretical and macro- 
scopic or statistical. The first scale, the microscopic correlation scale, characterizes 
the reactions taking place in the system, while the second scale measures the relax- 
ation of the system. For a weak coupling theory, in which we are interested in this 
paper, the former scale is much smaller than the latter scale.Q In a derivative expan- 
sion of some quantity F with respect to macroscopic spacetime coordinates X, we 
use the gradient approximation throughout: 

F(X, ...) ~ F(Y, ...) + (X - Yrd Y ,F(Y, ...) . (1) 

Let A(x,y) be a generic propagator. For the system of our concern, A(x,y), with 

x — y fixed, does not change appreciably in (x + y)/2. We refer to the first term on the 

2 It should be noted, however, that, as the system approaches the critical point of the phase 
transition, the microscopic correlation scale diverges. Thus, the formalism developed in this paper 
applies to the systems away from the critical point. 
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right-hand side (RHS) of Eq. (|1|) as the leading part (term) and to the second term as 
the gradient part (term). The self-energy part IT(a;, y) enjoys a similar property. Thus, 
as usual, we choose the relative coordinates x — y as the microscopic coordinates while 
X = (x + y)/2 as the macroscopic coordinates. A Fourier transform with respect to 
x — y yields 

A(x, y) = j^- A e- iP <^A(X; P) (X = (x + y)/2) (2) 

(P^ = (p°,p)) together with a similar formula for IT. The above observation shows 
that P M in Eq. (§) can be regarded as the momentum of the quasiparticle participat- 
ing in the microscopic reaction under consideration. We shall freely use A(x,y) or 
A(X;P) [U(x,y) or IT(X;P)], which we simply write A [II] whenever obvious from 
the context. 

The perturbative framework to be constructed accompanies the generalized Boltz- 
mann equation (GBE) for the number density of quasiparticles. The framework allows 
us to compute any reaction rate by using the reaction rate formula ||. Substituting 
the computed net production rates of quasiparticles into the GBE, one can deter- 
mine the number densities as functions of macroscopic spacetime coordinates, which 
describes the evolution of the system. 

For defmiteness, we take up QCD. However, the procedure has little to do with 
QCD and then the framework to be constructed may be used for any theory that 
includes gauge boson(s) with almost no modification. The plan of the paper is as 
follows: In Sec. II, on the basis of CTP formalism, we construct the bare gauge-boson 
and FP-ghost propagators, with which the "bare-X scheme" may be constructed. In 
Sec. Ill, we set up the basis for formulating the "physical-X scheme." In Sec. IV, 
we make up the self-energy-part resummed propagators. In Sec. V, we impose the 
condition that no large contribution appears in perturbative calculation, so that a 
"healthy" perturbative framework is constructed. It is shown that, on the energy- 
shell, the condition turns out to the generalized Boltzmann equation. In Sec. VI, 
we frame a concrete perturbative framework. Section VII is devoted to summary 
and outlook. Concrete derivation of various formula used in the text is made in 
Appendices. 
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2 Closed-time-path formalism 



2.1 Preliminary 

For definiteriess, we take up QCD and adopt the Coulomb gauge. We deal only with 
the gluon and the FP-ghost sectors in an out of equilibrium QCD system. Singling out 
the free parts of the gluon and of the FP-ghosts, we write the Lagrangian (density) 

as 

£{A {a) » (x) , r/ a) (x) , r/ (a) (x) , . . .) = - X - A^V^A (a)u + ?? (a) V^?/ a) + . . . , 

v „u = _ g ^Q2 + qvqv _ 1 UQP _ QO n H\ fff< _ Qp n »j\ . ( 3 ) 

a 1 J 
Here = (1,0), Aw? ("a" the color index) is the gluon fields, and 7/ ' and fjw 
are the FP-ghost fields. The CTP formalism is formulated || [7], by introducing 
an oriented closed-time path C (— C\ + C2) in a complex-time plane, that goes from 
— 00 to +00 (Ci) and then returns from +00 to —00 (C2). The real time formalism is 
achieved by doubling every degree of freedom: For example, for gauge fields, AW — > 
(4 a)M ,4 a)M ), where ^S^^cx) = A^(x ,x) with x e Ci and 4 a)/ >o,x) = 
iW^iOi x) with xo E CV A field with suffix 'z' (z = 1, 2) is called a type-z field. One 
can introduce a two-dimensional space, in which every field is a "vector" whose first 
(second) component is the type-1 (type-2) field. In equilibrium thermal field theory, 
this space is called the thermal space, so that we use this terminology in the following. 
The classical contour action that governs the dynamics of nonequilibrium systems is 
written in the form 



/ dx [ d^C{A ia ^{x),r ] ( - a) {x),f ] ia) {x),...) = [ ^dxo [ dx£(x) , 

JC J J~oo J 

£{x) = £{A^{x)^\x),^\x),..)-£{A^ (4) 



where stands for quark fields. £ here is sometimes called a hat-Lagrangian [cf. 
flOH . Throughout this paper, we do not deal with initial correlations (see, e.g., 0). 



2.1.1 Gluon sector 

Following standard procedure 0, the four kind of gluon propagators emerges: 

y)) = -iTv[T{A^{x)Af v {y)}p 
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{ab)iMi/ 
v 22 



— iTr 
— iTr 
— iTi 



A^(x)Ar(y)p 

T{A^(x)A?»(y)}p 



(5) 



where p is the density matrix, and T (T) is the time-ordering (antitime-ordering) sym- 
bol. We restrict ourselves to the case where p is color sing let, so that Ajf )(1U = 5 ab A%. 
Throughout in the sequel, we drop the color index. At the end of calculation we set 
Ai = A% |7[]. Let us use bold-face letter Aij(x,y) for denoting 4x4 matrix, whose 
{p z/)-component is A^(x,y). We further introduce a matrix propagator A(x,y), 
where the 'caret' denotes the 2x2 matrix in thermal space, whose (i j)-component 
is Aij(x,y). The matrix self-energy part Tl(x,y) is defined similarly. 

Deduction of the gluon propagator A is carried out below in subsection B. 



2.1.2 FP-ghost sector 

The FP-ghost fields are unphysical fields, so that they are absent in the system under 
consideration. Then, the FP-ghost propagators are the same as in vacuum theory. 
This means that the 2x2 matrix propagator that act on a thermal space is diagonal, 
whose diagonal elements are 



A g ab) (x,y) 
Af\x,y) 



ii 
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-i(0\T{ V { a \x)v?>\y)}\0) 

^(o\T{ v { 2 a \x)4 b \y)}\o) 



In what follows, as in the gluon sector, we drop the color indices. Setting rji = r]2 and 
Vi = V2, the FP-ghost propagator may be written, in 2 x 2 matrix notation, as 



A g (x,y)?3, 
d 4 P 
(2^ 



—% 



-iP-(x-y)_ 



p z 



(6) 



As will be shown in Appendix A, the FP-ghost self-energy-part matrix Il g is also 
diagonal, and enjoys a property [n g ] n + [n g ] 22 = with Im[n s (X; P)]u — 0. Then 
we may write 

U g (x,y) = U g (x,y)r 3 . 
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A f[ g -resummed FP-ghost propagator G g (= G g f 3 ) is defined through a Schwinger- 
Dyson equation: 



(■, \ • A„-ii (/ -r; (/ a„ • (;, r \\, r x r 

Here we have used the short-hand notation F ■ G, which is a function whose u (xy)- 
component" is 

[F-G](x,y) = Jd A zF(x,z)G(z,y). (7) 
Solving this to the gradient approximation, we obtain 

G g (X;P) = G g °\X;P) + G g 1 \X;P), 

9 1 ' ' 2 dP" dX^ [ 9 { ' h ' 

As will be observed in Sec. IV, to be consistent with the approximation we are taking, 
Gg 1 ' may be ignored, so that 

G 9 (X;P)r 



p 2 + n 9 (X;P) 
2.1.3 Vertex factors 

The vertex factors may be read off from Eq. (|4|): There is no vertex that mixes type-1 
fields with type-2 fields. Every vertex factor for interactions between the type-1 fields 
is the same as in vacuum theory, while the every vertex factor for the type-2 fields is 
of opposite sign to the corresponding vertex factor for the type-1 fields. 
In the next subsection, we deduce the gluon propagator, Eq. ([5]). 

2.2 Gluon propagator 

We compute Eq. (^) with A± = A 1 ^ (= A M ). For A M , we use the plane- wave decompo- 
sition with helicity basis in vacuum theory. Straightforward but lengthy calculation 
yields 

K^(x,y) = j ^e- iP <^k^(X-P) (X = (x + y)/2), (8) 
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where 
A^(X;P) 

Af (X;P) 
Af(X;P) 

Af (X;P) 



Af (X; P) + Af (X; P) + Af (X; P) + Af (X; P) + Af (X; P) , 

(9) 



-VP(p) A RA (P) + [A R (P) - A A (P)} f(X; P)M 4 



n tk n v - a—P^P u 

p2 



1 

p'2 



(10) 
(11) 



' (1+ =m(^v k) 



p 



2p(l-(^) 2 ) 

(l - ( P y) 2 ) V 



[A R (P)-A A (P)]/(X;P)M 4 



Af(X;P) 
Af(X;P) 



ig^g»j € vkpk N _(P) [Afl(P) - A A (P)] e(p°)M + , 

2p l p J 



(12) 
(13) 



(1 -2^)^ 



1 - (pf) 2 
^[p*(p x e 2 ^' +p j (p x e2 0*] 



ReX iR (P) 



1 - (j^) 2 

-(p x e^)^ - (p x e y ) j S iy ] e(p°)lmN LR (P)} 
x[A R (P)-A A (P)]e(p°)M + , 



(14) 



where the Greek letters, \x and u, take 0-3, while the Latin letters, % and j, take 1, 
2, and 3. To avoid possible confusion, we have written (x,y,z) for (1,2,3). In the 

1, and 



above equations, e y = (0, 1, 0), = d/dX\ e 123 = e xyz 



P 



M. 



± 



f(X;P) 
N(X;P) 



= -gffrif-Pff) (p = p/p), 
= p-(p. e ")e«, 

= Vk-(V X -P)p% 

Afl(P) 
A*(P)-A A (P) -A A (P) 

1 ±1 
±1 1 

e(po)X(X;P)-0(-p o ), 

i [X ++ (X; |p°|,e(p°)p) + X__(X; |p°|,e(p°)p) 
6(p°)n(X; \p°\,p) + 6(-p°)n(X- t \p°\, -p) , 



(15) 



A^)(P) 



1 



P 2 ± ip + ' 

(16) 
(17) 



(18) 



AL(X;P) 

N LR (P) = 
with 

N^(X; \Po\,i>] 



N ++ (X; \p°\,e(p°)p) - AL_(X; |p°|,e(/)p) 
AL + (|p |,e(p )p) 



2 

AL(P), 



(19) 
(20) 



d 3 qe~ 



■x 



Tr 



4(p-q/2)a c (p + q/2)p (£, C 



(21) 



Here |po| = Pi and Q = (q°, q) with g = q • p/f>- a>+ (a_) is an annihilation operator 
of a gluon with + (— ) helicity. a\_ (a^L) is a corresponding creation operator. In 
deriving Eqs. ([|) - ([2~TD, we have assumed that Tr[ap(p)a CT (q) p] ~ 0. We have also 
assumed that AL, Eq. (|19"D , and Nlr, Eq. (|^), are small quantities, so that A^ 



and A 3 U are of the order of gradient term Ai U that includes Vx^V- More precisely 
iV_ and N LR (P) are of order of (Vx-W) /p, and then, to the gradient approximation, 
AT-dependences of AL and of Nlr have been ignored. 

We observe that the structure of A^(P) in Eq. (|Tl|) is rather simple, which 
represents the propagator between "longitudinal gluons" as well as the gauge part. 
This is a reflection of the fact that, in the Coulomb gauge, only transverse modes are 
physical modes. Then, as in the FP-ghost propagator in Eq. (||), the (1 l)-component 
of A^, [A^(P)]n, is the same as in vacuum theory. 

From Eq. (|2l|) , follows 

9(±p°)P ■ d x N^(X; \p°\,±i>) = (|po| = P) ■ (22) 

Using this in Eqs. ([17]) and (|18|), we have 

P-d x f(X;P) = (\ Po \=p). (23) 

As is obvious from the construction, or as can be directly shown from Eq. ([23]), we 
see that 

- f 3 V^A pv (x, y) = -A»r(x, y)V pu f 3 = 6^5 4 (x - y) , (24) 

where T> w is as in Eq. ([|). Two equations in Eq. ( |2"2] ) are "free Boltzmann equa- 
tions." One can construct a perturbation theory in a similar manner as in ||, where 
a complex-scalar field system is treated. We call the perturbation theory thus con- 
structed the bare- A" scheme, since N obeys the "free Boltzmann equation." This 
scheme is equivalent [H to the physical-^ scheme, to which we now turn. 
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3 Construction of the physical- N scheme 



We aim to construct a scheme in terms of the number density that is as close as 
possible to the physical number density. To this end, first of all, we abandon the 
"free Boltzmann equation" fl23|). This means that / and A in the present (physical- N) 
scheme differs from the bare-iV scheme counterparts. Specification of / is postponed 
until Sec. VI. 

Now, in contrast to Eq. (24|), A pl/ is not an inverse of —r^D^. Straightforward 
calculation within the gradient approximation yields 

d A P 

x [A R (P) - A A (P)} PF(P) (P ■ d x f(X; P)) . 



r 3 V^ p A^(x,y) = -g^6 A (x-y)-ir 3 M + j^e- iP ^ 



In obtaining this, we have used P 2 [A R (P 2 ) - A A (P 2 )} oc P 2 5(P 2 ) = 0. Our procedure 
of constructing a consistent scheme is as follows: We further modify A p by adding a 
suitable A^ dd to A pv . The conditions for A^ dd to be satisfied are 

• A^ dd vanishes in the bare-iV scheme. 

• To the gradient approximation, the counterpart of Eq. ([24]) turns out to 

-r 3 XV - (X) J • (A pu + AZ d ) = 5;i , (25) 

where use has been made of the short-hand notation (]7|). '1' is the matrix in 
spacetime whose (x j/)-component is 6 4 (x — y), and L c is some (4 x 4) <g> (2 x 2) 
matrix function. 



It is a straightforward task to obtain the form of the required A 



in/ 

add - 



A^ d (X; .P) = iM + VF(p) [A 2 R (P) + A A (P)] P ■ d x f(X; P) , (26) 

from which we obtain for L pu in Eq. (|25| ). 

i^(x,y) = LfM_ 

= 2zM_| ^L e ~^-y) V ^ip)p . dxf{X -P) , 

with M_ as in Eq. ([16]). In obtaining Eq. (|25|) with A^ d as in Eq. (|26"D, we have 
used (Lc)^ ■ Aq P ~ {LcYv • (A up + A^ dd ), since the difference can be ignored to the 
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add J 



-r 3 l V 



gradient approximation. In a similar manner, we find that (A w + A 
— {L c )pu) = 8^1. Thus we have found that the 2x2 matrix propagator (A^ + A^ d ) 
is an inverse of (— f 3 V^ u l — (L C ) M „), so that the free action isQ 



(27) 



-- j d 4 xd 4 y%(x) [hV^5\x -y) + L^(x,y)M-\ A u {y) 



Since the term with L c (x,y) in Eq. (p7 ) is absent in the original action, we should 
introduce the counter action to compensate it, 



which yields a vertex 
iL^(x,y)M_ 



- j d A xd A y%(x)L^(x,y)M„A u (y) 



-2M_ J ^L e -* p <*-y^(p)P ■ d x f(X; P) 



(2* 



= ifl^(x,y). 



(29) 



For completeness, we construct a L c -resummed propagator. Since M + M_ = 
M_M + = and fl^ c ^ u oc only piece that changes the form under the L c - 

resummation is — Vj" {p)A RA in Ajf , Eq. ( [TD| ) with Eq. 



-pr(p)As um (^ 2/) = -^r(p) A RA (x, v ) - Y, [^ra {-n (c) • A^} 

I n=l 

= -VP(p) {A&(s,y) - [Ara ■ fl (c) • A RA ] (x,y)} 



A RA {x,y) + 2iM^ 



d 4 P 



(2tt)< 



-iP-(x-y) 



xA R (P)A A (P)P-d x f(X;P)} . 



Here, use has been made of A R / A M_A R ' A oc M + . Then, substituting (^A^ um ^ for 



A R U A , we obtain, with obvious notation, 



(A~(X;P))^ + A^ d (X;P) 

= -VF(j>) { [A RA (X; P) + M + [A R (P) - A A (P)} f(X; P) 
-iM + VF(p) [A R (P) - A A (P)} 2 P ■ d x f(X; P)} . 



(30) 



3 Note that the Lagrangian density corresponding to the term with L£ u {x, y) in Eq. (g7|) is nonlocal 
not only in 'space' but also in 'time.' Here it is worth mentioning the so-called |po (-prescription. 
With this prescription, at an intermediate stage, we have (L' c ) tJ ' 1 ' {x,y), which is local in time. [See 
§, 1 for details.! 
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In closing this section, we emphasize that f(X; P) in the present scheme is an ar- 
bitrary function, provided that /(X° n ,X;P) [= -0(-p°) + e(jP)N(X? n ,X;P)], with 
X™ the initial time, is a given initial data. We have introduced the counteraction 
A c , Eq. fl28D , so as to remain on the original theory. Thus, it cannot be overempha- 
sized that the schemes with different /'s are mutually equivalent. If we choose the / 
(= fs) that subjects to the "free Boltzmann equation" fl23p, the scheme reduces to 
the bare-iV scheme in the last section. In Sec. VI, we shall choose /, with which a 
well-defined perturbation theory is formulated. In any case, it is natural to choose /, 
so that, when interactions are switched off, / turns out to 



4 Resummation of the self-energy part 
4.1 Preliminary 

Throughout in the sequel, we restrict to the strict Coulomb gauge a = [cf. Eq. (|3|)]. 
The bare propagator consists of six pieces, At + Al + 2f =1 Aj + A a dd, Eqs. (|) 
and (HH) . At + Ax is the leading part and Ai and A a dd are the gradient parts, 
which represents variation in the macroscopic spacetime coordinates X^, through 
first-order derivative dx^f{X; P). We have assumed that A 2 and A 3 are of the same 
order of magnitude as the gradient parts. Interactions among the fields give rise to 
reactions taking place in a system, which, in turn, causes a nontrivial change in the 
number density of quasiparticles. Thus, the self-energy part n ties with the gradient 
parts. More precisely, n is of the same order of magnitude as A ad d (A t + A L ) 

and Aj (At + A^j (i = 1—3). Hence, in computing n in the approximation 
under consideration, it is sufficient to keep the leading part (i.e., the part with no 
X M -derivative). Then tl(X; P) may be decomposed as 

W V {X;P) = V T v (p)fl T (X;P)+n^n v fl L (X;P) 

+ El (jfpL + pa n -) n c (X; P) + pH^-fl D (X; P) . (31) 
P 

Here and in the following, we use such a notation as for denoting a four-vector 
whose Oth component is zero: = (0, v). Then, = (0, p). 

Within the gradient approximation, it is sufficient to perform a II-resummation 
for the leading part, At + At- This is because the corrections to other pieces due to 
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the resummation are of higher order. Thus, for the gradient parts of A as well as for 
A a( jd, one can use the formulae in the bare-iV scheme in Sec. II [cf. above argument 
after Eq. ([30])]. In particular, for / in the gradient parts, one can use / (= /#) as 
in Eq. ( p3|) in the bare-iV scheme. For the purpose of performing resummation, as 
usual, it is convenient to introduce the "standard form" [n| [3], f|, |!| 



Bl ■ Adiag • B 



diag 



d z x Id z 2 B L (x, zx)A d i a , s (z 1 ,z 2 )B R (z2,y) , (32) 



Bi 




pP pv 




P v n^ 



p2 
,2 



p2 



(33) 



where / (= f(x, y)) is the inverse Fourier transform of /(X; P). It is to be noted that 



-A R u , A Jx — y) is an inverse of Zy, Eq. 



U P^R(A) 



AMP Jjv _ 
R(A) U p~ 



pv 



(34) 



Computing Eq. ( p2[ ) to the gradient approximation, we obtain 



A 



pv 
base 



%y) = &r(x,y) + A^(x,y) + A^ d (x,y) 



b l ■ A^ ag ■ B R 



{x,y) + M+A%(x,y) 



A%(X;P) 



p V 



^(^Vi + ^) / + 2^p.V x / 



p- 



(35) 
(36) 



where P/P 2 is the principal part of 1/{P 2 ± i0 + ). As mentioned above, within the 
gradient approximation, it is sufficient to take At + A^ as a "resummed part." It 
is obvious from the above observation that one can freely include gradient part(s) 
into the "resummed part." We include the gradient part A^ d and take A^ se as the 
"resummed part." 

It is to be noted that, from Eqs. (0), (|TT|) , ([IB]) , and fl2"6|), follows 



EH' i+J (A b asck = 0. 



(37) 



This relation holds for A + A a dd (in place of Ab ase )- 
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4.2 Self-energy-part resummed propagator 

As will be shown in Appendix A [cf. Eq. (|A.7|) 1, 



E n *i = ( 38 ) 
holds. A II-resummed propagator G obeys the Schwinger- Dyson equation: 

G = Abase — Abase • II • G = Abase — G ■ II ■ Abase • (39) 

It is clear that, for solving this equation, it is convenient to introduce the U B- 
transformed" quantities 

G = B L l -G.B- R \ (40) 

A base = B L l ■ A basc • B}? = I Ar A * I , (41) 

V -A A 



U = B R -U-B L 



o -u A 



Here A^ is as in Eq. (|36D , and II_r(a) and II ^ are, in respective order, the inverse 
Fourier transforms of 

U R{A) {X;P) = n n + n 12(21) (42) 
U K (X;P) = [l + f(X;P)]U 12 (X;P)-f(X;P)U 21 (X;P) 

-1{£ R + £ A , /}pg . (43) 

In Eq. fl43|), {..., ...}rb. is defined as 

i a m - dA ^ p ) dB ( x > p ) dA ( x > p ) dB ( x > p ) (AA\ 
1 ' lpB - " dx» op, dP^ dx» • 1 ] 

In obtaining Eq. ([y]) [Eqs. ( fffi) and (f43f )1, use has been made of Eq. (|37D [Eq. (|38|)1- 
Although the last term in Eq. fl43|) may be dropped to the approximation under 
consideration, we have kept it. 

Now Eq. ( |3"9"D is transformed into 

Q- .A^base .A.base ' U ' Gl .A_b ase Gr • II • .A_b ase , 
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which can easily be solved to give 



G 

Gr(a) 
G K 



G R 


^R(A) 



Gk 
-G A 



+ n 



R(A) 



G R ■ [Al 1 ■ A K ■ A^ 1 + n^J • G A . 

Substituting this back into Eq. ([40]), we obtain, after Fourier transformation, 
G(X;P) ~ G RA (X-P) + M + [G R (X;P)-G A (X;P)}f(X;P) 



(45) 
(46) 



+M 4 



G K (X;P)--{G R + G A , /}pg 



where 



Gr A (X; P) 



G r (X;P) 
G r {X-P)-G a {X-P) 

Here we note that II consists of two pieces, 

ri = ri loop + ri (c) , 



o 

-G A (X;P) 



(47) 



(4$ 



(49) 



where n loop is the contribution from loop diagrams and n( c ) is as in Eq. (|29"D , which 
has come from the counter action A c . It should be remarked that some n loop contains 
internal vertex(es) itl^ c \ [cf. the argument at Sec VI below.] From the above 
definitions of T1 R ( A ), Tlx, and & c \ we have 



n 



R 



n 



o 



n 



K 



(X;P) 



-2iV^{i>)P-d x f{X;P). 



(50) 



From Eq. (|45|), we obtain, after some manipulation, 



G R{A) {X-P) ~ Gf {A) {X-P) + G^ [A) {X-P) 

R(A) \ i / - 



T IP J 



/j, v 



+ 



P 2 -n T R{A) (X;P) p> + IL& A) {X;P)' 



(51) 



R(A)\ 



G R(A)( X ' P ) 



2p 2 p2 _ U T 



R(A) 



P2 _ r\T 



if*\ - iM) n 



fl(A) 



Po 



p 2 + n L 



(52) 
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Although G { ^ {A) can be ignored to the gradient approximation, we have displayed it. 



Straightforward calculation of the last term in Eq. ( f4T| ) using the definitions (|45|), 
(M), and (M) yields (cf. Appendix B) 



where 



G 



K 



G 



(2) 
K 



where 



G K (X;P)--{G R + G A , /} pa 



n (l) n (2) 



— I 



—% 



t yoop 



2P.a x / + {Ren^, /} pB +i(n 
2p-v x / + {Ren^, /} pB +z(n^) l0 °p 



T 



(P 2 -nS(P 2 -n^ 



Po 



-ImlTg 



(p 2 + n^)( P 2 + n^ 



,(53) 



(P 2 -iW + rri) (P 2 -iW + n£) 
^(p^nl---)(^i/-^/) 

(P»- r ng). ("'•*/ + {"£■ /} P J 



-MRe 



(54) 



(n^) loop .[i + /](n^) loop -/(n^ 

In obtaining Eqs. (^) and (|54"D, use has been made of (ll 
proved in Appendix A. 



loop 



loop 



T(L) 
A 



(55) 

11^ , which is 



5 Generalized Boltzmann equation 

5.1 Energy-shell and physical number densities 

For later use, referring to Eq. fl5"T|), we define the energy-shell for the transverse mode 
by 



Re 



P 2 -Ren£(X;P) 



0. 



(56) 



p°=±u> T (X;±p) 

It is well known |2[] that, in equilibrium quark-gluon plasma, "longitudinal" mode 
called plasmon appears for soft p (= O(gT)) [g the QCD coupling constant]. The 
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energy-shell of such modes is defined by [cf. Eq. 



Re 



p 2 + Ren£ (X; P) 



0. 



(57) 



JpO=±^ L (X;=Fp) 

Useful formulae that hold on the energy-shell are displayed in Appendix C. 

We restrict to the po > part. The po < part yields the same result, since a 
gluon and a corresponding antigluon is the same. In some of the formulae in this 
subsection, we shall drop the argument X. 

For obtaining expressions for the number densities of the transverse and "longi- 
tudinal" modes, we compute the statistical average of the operators, 

d 



Mr(x) = -iAf\x) — A { r\ x ) 

J OXn J 



M L {x) 



d 
dx 



A { -\x) 



Here the superscript "(+)" ["(—)"] stands for the positive [negative] frequency part. 
Statistical averages of At and of Ml yield, in respective order, 

d 4 P 



Tr [M T (x)p] 
Tr [N L {x)p] 



(27T) 4 

d A P 



—i 



(27T) 



0(Po)Po [G{i(X;P) + Gi{(X;P) 
0( Po )Po \G° l l(X;P) + G 2 1 (X;P) 



(5? 



It is to be noted that "x" here is a macroscopic spacetime coordinates [cf. Sec. I]. 

We first compute the leading contributions to Eq. (|58[ ). Substituting the leading 
parts of G 2 i and of G\ 2 (cf. Eq. (pETD) and using Eq. fl5"T|), we obtain 



Tr [M T {x)p\ = 2i[ 
Tr[Af L (x)p] = - 



d A P 

d 4 p 



%o)Po 
6{Po)Po 



1 



Pi-Til 



c.c. 



— c.c. 



The narrow-width approximation^ ImiT^ 



T(L) 



P 2 + n£ 

-> — e(po)0 + yields 



[2/(X;P) + l] 
[2/(X;P) + l] 



Tr [N T {x)p\ 
Tr [U L {x)p] 



d 3 p 
(2tt) 3 
d 3 p 



Z T (w T (p),p)n(a; T (p),p) 



f d p 



(59) 
(60) 



In the case of equilibrium system, the narrow-width approximation is a good approximation for 
hard modes but not for soft modes. 
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Here n is as in Eq. (|Hf ), and stands for the contribution from 2f + 1 3 e(p°) [ cf. 
Eq. flI7p], which is the vacuum-theory contribution corrected by the medium effect. 
Z's in Eqs. (^) and (|60|) are the wave-function renormalization factors, Eqs. (|C.1| ) 
and ( |C.4j ) in Appendix C. If there are several modes, summation should be taken over 
all modes. The factor '2' on the RHS of Eq. (59) is a reflection of the two independent 
polarizations. 

Thus, we have learned that n is the number density of gluonic quasiparticle. 
Undoing the narrow-width approximation yields further corrections to the number 
densities. 

Let us turn to analyze the contributions from the other parts of G in Eq. (|47j). 
Inspection of Eq. fl5|) with Eqs. (g7[) - shows that all but G K \ Eq. (0), yield 
well-defined corrections to the physical number densities due to the medium effect. 
G K contains 

(61) 



[j» - n@ [p 2 - (n^ 

and 

(62) 



In the narrow-width approximation ImlT^ — ► — e(po)0 + [ImlT^ — > —e(po)0 + ], Eq. (plj) 
[Eq. (|62"D 1 develops pinch singularity at p = ujt(p) [po = ^l^p)] in a complex p°- 
plane. Then the contributions of Gx to Eq. (|58|) diverge in this approximation. 
In practice, lmIl^ L ^ (oc (? 2 ) is a small quantity (at least for hard modes), so that 
the contribution, although not divergent, is large. This invalidates the perturbative 
scheme and a sort of "renormalization" is necessary for the number densities [^, [|, [5] . 
This observation leads us to introduce the condition G K = on the energy-shells: 

2P . Bx} + {Ren- /} pB + i(^] K ^ X:±p} = , (63) 
2p • V X / + {Reni /} pB + Hn^] K ^ L(X;±p) = , (64) 

These serve as determining equations for so far arbitrary /. (See below, for more de- 
tails.) Now the above-mentioned large contributions, which turn out to the diverging 
contributions (due to the pinch singularities) in the narrow-width approximation, do 
not appear. Thus, the contributions from G K to Eq. (p8|) also yields well-defined 
corrections to the number densities. 
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5.2 Generalized Boltzmann equation 



We are now in a position to disclose the physical meaning of Eqs. fl63|) and (|64]) . With 
the help of the formulae in Appendix C, Eq. (|63D with p — > u>t(X; p) turns out to 



(Z T (X;P)) 



+ 



-i 



9 / v \ 9 

^ + v T (X;p).— 



n(X;P) 



<9RelT5(X; P) dn(X; P) 



2u; T (X; p) dX^ 



(n 



loop 



dP» 
(X;P) 



(65) 



2uj T (X;p) 

Here (= dur^X; p)/dp) is the group velocity of the mode [cf. Eq. 00.20 1. As will 
be shown in Appendix D, the RHS of Eq. ( |65"D is related to the net production rate, 
-^net P j °f the mode p° = Ut(X;p). Using Eqs. ( p.l| ) and ( |0.3|) , we obtain 



9 , v ^ 9 

^ + vt(X; P ).— 



n(X;P) 
duj T (X;p)dn(X;P) 



dX» 

This can further be rewritten in the form, 



dP„ 



net p 



(A';p) 



Similarly, Eq. 



duj>i>{X; p) dn 
dX dp = 

with p = Ul(X; p) yields 

(^ + v,(X;p).A 

8ujl{X; p) dn 
dX dp ~ 



net p 



(X;p). 



(66) 



n(X;^ L (X;p),p) 



rL P (X; P ) 



(67) 



n here is essentially (the main part of) the relativistic Wigner function, and Eqs. ( |66"D 
and ( pTQ are the generalized relativistic Boltzmann equation for gluonic quasiparticles 
[cf. 0]. 

For the sake of defmiteness, we are taking QCD and have dealt with the gluon 
sector. The quark sector has already been dealt with in ||. As far as the gluon 
sector is concerned to the gradient approximation, one can use the "free quark" 
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distribution function that subjects to the "free Boltzmann equation." This can be 
seen as follows: In the gluon sector, the quark distribution function enters through 
the quark propagators involved in the gluon self-energy part II. As observed at the 
beginning of Sec. IV, II is related to the gradient part of the gluon propagator. Thus, 
in calculating II, one can use the Oth-order expression (of the derivative expansion) 
for the propagators involved in n. 



6 Perturbation theory 

As has been discussed in the preceding section, the propagator in the physical-iV 
scheme is free from the pinch singular term (in the narrow-width approximation) 
and then the perturbative calculation of some quantity yields "healthy" perturba- 
tive corrections. For constructing a concrete perturbative scheme, one more step is 
necessary. 

To extend the conditions ( |63D and (|64D to off the energy-shells, we divide / into 
two pieces [cf. Eq. (|17|)1 

f(X;P) = e(p°)N(X;P)-6(-p ) 
= f(X;P) + f (X;P) 

= [e(p )N(X;P)-9(-p ))+e(p )N (X;P). (68) 

fo (and then also /) is defined as follows: Let lZi(X;p) (i = 1,2,...) be a region in 
a p°-plane that includes zth energy-shell. We choose 7Zi(X; p), such that, for i ^ j, 
IZi fl IZj = 0. On each energy-shell, N (X;P) = N(X;P), and, in whole p°-region 
but Ki(X;p) (i = 1,2,...), N (X;P) vanishes. dN (X;P)/dX and dN (X;P)/dP 
exist and iV (X; P) obeys 

2P ■ d x N + {ReU T R , N } pB = -i , 

2p • VxiVo + {ReIT£, N } pB = -i (n£-) l0 ° P . (69) 

Then, G^^ u in Eq. (B3I) turns out to 



G K 



2P.^/ + {Ren^, /} +z(n 



T yoop 



V; 



(P 2_nT) ( p2_n 



2p • V X / + {Ren^, /} pa + » ( n' )>-] (p2 + uLR)(p2 + JiLA) .(70) 
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It is obvious from the above construction that this G^'^ u does not possess pinch 
singularities in narrow-width approximation, and thus "healthy" perturbation theory 
is established. 

It is worth mentioning that there is arbitrariness in the choice of the regions 
TZi(X;p) (i = 1,2,...). Furthermore, the choice of the functional forms of / and of 
fo is also arbitrary, provided that 

f(X° = Xl, X; P) + f (X° = Xl, X; P) = f(X° = X° n , X; P) , 

where f(X° = X® n ,X.;P) is the initial data with X® n the initial time. As has been 
discussed at the end of Sec. Ill, these arbitrariness are not the matter. 

To summarize, to the gradient approximation, the (resummed) propagator Q^ v of 
the theory is 

= & v + Af + Af + Af . (71) 

Here Af - Af are as in Eqs. QT2|) - (|TJ) and & v is as in Eq. (g7p with Eqs. (gSJ) - 
(|55l), provided that G^^ is given by Eq. ([ft]). / consists of two pieces as in Eq. (|68|) . 
fo (— e (p°)No) subjects to Eq. (|69|) , which is to be solved under a given initial data. 
It is to be noted that Af - A3 v are the gradient parts, so that, if one wants, for 
/ in Af - Af , one can substitute the solution to the "free Boltzmann equation," 
Eq. (£3). 

Determination of / or N proceeds order by order in perturbation theory to get 
/ = fo + fi + f2 + •••• Provided that perturbative computation is completed up to 
(n — l)th order, one can proceed to nth order calculation, f or N enters through 
Af se , Eq. and through tl^ (oc P ■ df), Eq. (^). For / in Af se and in 

jj( c ) j n Eq, (^), we substitute 2" =0 /j, while for / in FI^ involved in some FI loop , 
Eq. (|49|), we use appropriate Yh=o fi (j ^ n ~ !)■ Then, proceeding as above, one can 
determine f n , and thereby any physical quantity may be computed. Here we recall 
that we are taking the gradient approximation. Then, in proceeding to higher orders, 
one should check whether or not the consequences of the calculation of the order under 
consideration is consistent with the gradient approximation. One suspects that, in 
general, only first few orders meet this criterion. 

As to the vertex, as has been mentioned at Sec. II A, the matrix vertex V is 
diagonal and its (2 2)-component is of opposite sign relative to the (1 l)-component. 
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It is to be noted that the two-point vertex iL^ u (x, y)M- in Eq. (p9|) has already been 
built into G^ u (cf. Eqs. ( pEDD and floT)|)) and is absent in the perturbative framework 
using & v in Eq. (0). 

7 Summary and outlook 

In this paper we have dealt with out-of-equilibrium perturbation theory for gauge 
bosons and FP-ghosts in Coulomb gauge. The gauge-boson and FP-ghost propagators 
are constructed from first principles. Only approximation we have employed is the 
so-called gradient approximation, so that the perturbative framework applies to the 
quasiuniform systems near equilibrium or the nonequilibrium quasistationary systems. 
The framework allows us to compute any reaction rates. 

There comes out naturally the generalized Boltzmann equation (GBE) that de- 
scribes the spacetime evolution of the number densities of quasiparticles, through 
which the evolution of the system is described. 

The GBE for a quark-gluon plasma (nonequilibrium QCD) is "directly" derived in 
|T2|1 [cf., also, Ul3| ]. The "derivation" of the GBE in this paper is quite different from 
that in |TJ] . What we have shown here is that the requirement of the absence of large 
contributions from the perturbative framework leads to the GBE. This means that 
the quasiparticles thus defined are the well-defined modes in the medium. Conversely, 
if we start with defining the quasiparticles such that their number density functions 
subject to the GBE, then, on the basis of them, well-defined perturbation theory may 
be constructed. 
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Appendix A: Properties of the self-energy parts 



A.l: Gluon self-energy part 



Let us start with studying the properties of the propagator. From Eqs. 
Eq. (p6|) , we obtain 



SOf) and 



iA^ (21) (X;P) 



iA%(X;P) = iA%(X;-P) 



^12(21) 



(X;P)=A^ (12) (X;-P) 



where use has been made of f(X; P) = —9(—p ) + e(po)n(X; e(p )p). 

Now, we observe the properties of vertices. As has been mentioned in Sec. IIA 
and Sec. VI, the vertex matrix is diagonal and its (2 2)-component is of opposite sign 



relative to the (1 l)-component, V = diag(t> 



(1) w (2) 



diag(«,-v). For QCD, let 



V3 = (vs(P), —vs(P)) be a 3-gluon-vertex factor and V4 = (^4,-^4) be a 4-gluon- 
vertex factor. Noting that v^(P) is real and linear in P, and t>4, being independent 
of P, is pure imaginary, we see that 



v£\P) 



v 3 (P) = -v 3 (-P)=v^(-P) 



(2) 

—V 4 = V\ 



holds. There is one more (two-point) vertex (|29|), which satisfies 



iU?{X\ P)M- = iL^(X; -P)M- . 



(a.i; 



(A.2) 



For quark |J and FP-ghost, one can write down the similar relations, which we do 
not display here. 

From the diagrammatic analysis of tl(X; P) using the above relations, one can 
straightforwardly obtain the relations, 



^12(21) » 



[n^(X;P)}* = -Hg(X;-P) = -Il%{X;P), 

,(X;P)]* = -n^ (12) (X;-P) = -n^ (21) (X;P). 

Using the decomposition (0H) for the leading parts of II, we obtain 

[(n s (x ; p)) n ]* = -(n 5 (X;P)) 22 , 



(Ks(X;P)) 



12(21) 



12(21) ' 

(S = T,L,C,D) 



(A.3) 
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The last equation shows that (ILs(X; -P)) 12 (2i) are pure imaginary. 
Let us turn to analyze II in configuration space: 



II(x,y) 



Y[ [dzi dui dvi 



i=l,2„ 



x n(z, y; zi, z 2 , ...; (u u (tt 2 , ua), •••) 



(A.4) 



where x and y are the spacetime points to which the external legs are attached, z% 
ii = 1,2,...) are the internal vertex-points, and (v,i,Vi) [i = 1,2,...] are the sets of 
internal vertex-points of iL c (ui,Vi)M_. From Eq. @, which leads to Eqs. (§) - (|20|) , 
and Eq. (|^), we see that 



A u (x,y) = A 2i (x,y) for x° > y° 



1,2). 



(A.5) 



Here we use the "largest-time technique" |]14|| . Using the relation ( |A.5| ) and Eq. ( |A.1| ) 
(in configuration space) together with their quark and FP-ghost counterparts, one can 
show that, in Eq. (|A.4| ), when z\ of z®, ...) is the largest time, then tl(x,y) = 0. 
As can be shown by using Eq. ( |A.2|) in configuration space, this is also the case for 
the largest w° or for the largest v® of v®), (u®, ^2), ■■•]• Then, in Eq. ( |A.4j ), the 
largest time is x° or y°. 

Similar argument as above shows that 



n 2 i(x,y) = -Il u (x,y), n 22 (x,y) = -n 12 (x,y) 
Un(x,y) = -U 12 (x,y), U 21 (x, y) = -U 22 (x, y) 



for x° > y° , 
for x° <y°. (A.6) 



From the first (second) relation, we get 11^(2, y) [= IIn + n 2 i] = for a; > y°, and 
Hr(x, y) [= IIn + II12] = for x° < y°, as they should be. From Eq. (|A.6|) , follows 



U ij (x,y) = 0. 

Using the decomposition (pT|), after Fourier transformation, we obtain 

j2 (n 5 (X;P)) ij= (S = T,L,C,D). 

The relation ( |A.3| ) leads to 

U S A (X; P)] * = U S R (X; P) (S = T, L, C, D) . 



(A.7) 



(A.f 



(A.9) 
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A. 2: FP-ghost self-energy part 

In every diagram representing the FP-ghost self-energy part fl g , two external FP- 
ghost lines are connected with each other. Then, the fact that the bare FP-ghost 
propagator is diagonal [cf. Eq. (H)] tells us that the matrix self-energy part Yl g is also 
diagonal, fl g = diag(n^, n< 2 >). 

In a similar manner to the above subsection, we obtain [cf. Eqs ( |A.8| ) and (|A.9|) 1 



n( 1 )(X;P) + nf(x ; p) = o, 

[n«(X;P)]* = -nf (X;P). 



Then, we see that 

fl 9 (X; P) = f 3 U g (X; P) , lmU g (X; P) = . 



Appendix B: On the derivation of in Eqs. ( |53f) 
and (gjj) 



Here we deal with a part of Gk in Eq. 



G R ■ A R ■ A K ■ A A ■ G A 



d 4 p 



e -P(x-y)G K (X-P) 



G K (X- P) = G R (X; P)A R L (P)A K (X; P)A A l (P)G A (X; P) , 



(b.i; 



which is valid within the approximation under consideration. We are adopting the 
strict Coulomb gauge a = 0. Since (A R \P)Y V = (A^(P))^ = -V U {P) [cf. 
Eq. (|34|)1 contains a term —g^ t g v ^p l p' /a (a — > 0), care should be taken for computing 
Eq. QB.lp . We need 0(a) contribution to G^). For the present purpose, it is 
sufficient to obtain the leading-order expression, which is obtained from Eq. (|45|): 



G R \ A) (X-P) ~ Gf^{X-P) + aG R r A) {X-P) 



P 



'p 2 + U c 



R(A) 



P \p + 1 - i -R(A) IP \ P + ii /?(.4) 
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where G^,^( is as in Eq. (|51] ) and p± = (0, p) as in Eq. (|31|). 
We observe that 



pi/ 



0' up n - — 



p 2 + Il c 



p 2 + n L 



R(A) v 

p Q n 

R(A) 



Substituting these and Eqs. ( pT|) and fl36|) into Eq. ( |B.1| ), we obtain, after some 
manipulations, 



G K {X-P) 



2p 2 n ll n v 



(p 2 + U L R )(p 2 + U L A 



-P-Vx/ + 



p 2 -it£ 



+ ™ ^ +P0 



> 2 + n£)(P 2 -n^) (p2 + n ^(p2-nS) 



Appendix C: On the energy shell 



Here we display some formulae, which hold on the energy-shells of quasiparticles [cf. 
Eqs. (|56| ) and (E7I) with Eq. ([51])]. In most formulae in this Appendix, the argument 
X is dropped. We restrict to the positive po part. 
Transverse modes: 

We first define the wave-function renormalization factors through taking derivative 
of Eq. fl56|) with respect to p°: 

1 <9ReIl£G» ,p) 



(Z T (u; T (p),p)) 1 = 1 - 



p° = u t (p) 



The group velocities of the modes are obtained from the definition (|56|): 

du*r(p) Z T (a; r (p),p) 



vtIPJ 



x 



dp wr(p) 
■ ldRelT£(p ,p) 



<9p 



P° = <^t(p) 



Differentiation of Eq. QSTf ) with respect to X leads to 

dur(X; p) Z T (u T (p),p) d Re IT^X; ^ r (X; p), p) 



OX 



dX 



(CI) 



(C.2) 



(C.3) 
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''Longitudinal" modes: 

From Eq. fl5T|), we obtain 

[Z L (u L (p),p)) 



-1 



1 dReU^(p°,p 



2cu L dp 



(C.4) 



v MPj 



duj L (p) 



dp 

Z L (u L (p),p) 
^l(p) 



P + 



l<9Ren£(p°,p 
2 dp 



P D = «i(p) 



Appendix D: Net production rates 

From Eqs. (f)TD, Q56|) and (|57|) , we see that the projection operators ontop = ojt (p° = 
c^x) mode is V^ u {p) (n^n u ). Then, the production (decay) rates of the transverse 
modes and of the "longitudinal" modes are written as |7|, 0, [|. |J 

i 1 

T p(d)(p) = 2uT ^ z TMp),p)-n^ i21) (v T (p))^ 



and 



r^)( P ) = ^%(p),p)iT 
" 27J7 L 12{21) ' 



12(21) n M n f 



-'L 

respectively. Here Z's are the wave-function renormalization factor, Eqs. ( |C.1| ) and 
T^j . Thus, the net production rate is 



SI(P) = [l + ^ T(L) (p),p)]r^)(X;p)-n(u; T(L) (p),p)rJ (L) (p) 

Z T(L) (cu T{L )(p),p) (n^ (L) ( WT(L) (X;p),p)) l0 ° P , (D.l) 



2ut(l) 

loop 



where (ll^} 1 ^ is as in Eq. fl55|). 



26 



References 



[1] J. D. Bjorken, Acta Phys. Polon. B 28, 2773 (1997); See also, e.g., K. Rajagopal, 
in Quark-Gluon Plasma 2, edited by R. C. Hwa (World Scientific, Singapore, 
1995), p.484. 

[2] M. Le Bellac, "Thermal Field Theory," Cambridge University Press, Cambridge, 
England, 1996. 

[3] A. Niegawa, Prog. Theor. Phys. 102, 1 (1999). 

[4] A. Niegawa, Phys. Rev. D 62, 125004 (2000). 

[5] A. Niegawa, [hep-ph/01012T6 . 



[6] J. Schwinger, J. Math. Phys. 2, 407 (1961); L. V. Keldysh, Zh. Eksp. Teor. Fiz. 
47, 1515 (1964) (Sov. Phys. JETP 20, 1018 (1965)); R. A. Craig, J. Math. Phys. 
9, 605 (1968); J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

[7] K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

[8] N. P. Landsman and Ch. G. van Weert, Phys. Rep. 145, 141 (1987). 

[9] A. Niegawa, K. Okano, and H. Ozaki, Phys. Rev. D 61, 056004 (2000). For the 
case of equilibrium systems, see, A. Niegawa, ibid. 57, 1379 (1998). 

[10] H. Umezawa, "Advanced Field Theory — Micro, Macro, and Thermal Physics." 
AIP, New York, 1993; Y. Yamanaka, H. Umezawa, K. Nakamura, and T. Arim- 
itsu, Int. J. Mod. Phys. A9, 1153 (1994); Y. Yamanaka and K. Nakamura, Mod. 
Phys. Lett. A9, 2879 (1994); H. Chu and H. Umezawa, Int. J. Mod. Phys. A9, 
1703 and 2363 (1994). 

[11] P. Zhuang and U. Heinz, Phys. Rev. D 57, 6525 (1998); S. Ochs and U. Heinz, 
Ann. Phys. (N.Y.) 266, 351 (1998). 

[12] J.-P. Blaizot and E. Iancu, Nucl. Phys. B557, 183 (1999). 

[13] S. Mrowczyhski, Phys. Part. Nucl. 30, 419 (1999) [|hep-ph/9805435|| , and earlier 
works quoted therein. 



27 



[14] R. Kobes, Phys. Rev. D 42, 562 (1990). 



28 



